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Abstract 

In the last few years, interest in radio detection of cosmic ray air showers has risen continuously. By now, large- 
scale application of the radio technique is under investigation in the framework of LOFAR and the Pierre Auger 
Observatory. The experimental efforts are accompanied by new approaches to describe and model the underlying 
radiation mechanisms, to lay the foundation for an interpretation of the experimental data. In this article, I review 
the current radio emission theory and simulations, with slight focus on the geosynchrotron model and its predictions 
for the information content of the radio signals. 

Key words: cosmic rays, extensive air showers, electromagnetic radiation from moving charges, computer modeling and 
simulation 

PACS: 96.50.S-, 96.50.sd, 41.60.-m, 07.05. Tp 



1. Introduction 



2. The early days 



Radio detection of cosmic ray air showers 
promises to nicely complement the existing particle 
and fluorescence detection techniques for extensive 
air showers (EAS). It offers 100% duty cycle, a 
calorimetric measurement, information about the 
longitudinal evolution of the shower and very good 
angular resolution. To exploit these advantages, 
however, the radio emission physics has to be un- 
derstood in great detail. Efforts to interpret radio 
signals from extensive air showers reach back to the 
early days of radio detection of cosmic rays in the 
mid-1960s. With the revival of the field in the last 
few years, interest in simulations and theory of radio 
emission from EAS has also been increasing contin- 
uously. After a very short mention of the historical 
approaches, I will give an overview of the theory 
and simulation activities of these last few years and 
the latest predictions on how to use radio detection 
to derive air shower parameters of interest. 
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In the year 1962, Askaryan for the first time sug- 
gested that charge excess in extensive air showers 
should produce coherent radio Cherenkov emission 
in the metre wavelength range. These theoretical 
predictions 0; 0] led to the first experiments and 
the successful detection of radio emission from EAS 
by Jelley et al. Q . After the successful detection, a 
flurry of activity ensued, both in experimental de- 
tection and theoretical studies of the emission mech- 
anism. In 1966, Kahn & Lerche Q proposed that 
the emission could also be of geomagnetic origin. 
Their model predicted the emission to arise from the 
transverse currents of electrons and positrons being 
deflected in the Earth's magnetic fleld and the asso- 
ciated dipole propagating through the atmosphere. 
Many more articles developing models of increasing 
sophistication were published in the following years 
(for a review please see Q), but the two main lines 
of thought remained those of Cherenkov radiation 
and geomagnetic radiation. When the field of radio 
detection of cosmic rays virtually ceased to exist in 
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the niid-1970s, the consensus, mostly based on po- 
larisation measurements of the emission, was that 
the radio emission was dominated by a geomagnetic 
effect. 

3. Modern radio theory 

In the year 2001, the experimental efforts for radio 
detection of cosmic rays with modern digital tech- 
nology began, and it was clear from the beginning 
that there was a need for more detailed theory and 
simulations of radio emission from EAS. Based on 
the historical results, the expectation was that the 
radio emission was dominated by a geomagnetic ef- 
fect. 

3.1. Frequency- domain analytical geosynchrotron 
model 

A new approach at the problem of radio emis- 
sion from cosmic ray air showers was taken with the 
geosynchrotron model first proposed in Q and then 
worked out in detail by Huege & Falcke in a num- 
ber of publications 0; B 0] ■ In addition, the geosyn- 
chrotron approach was explored by Suprun, Gorham 
and Rosner [l3|, who however did not follow the 
topic further. The geosynchrotron model describes 
the radio emission as coherent, synchrotron-like ra- 
diation from the electrons and positrons being de- 
flected in the Earth's magnetic field. In a first step, 
the emission was calculated with an analytical ap- 
proach in the frequency domain, based on parame- 
terisations of the relevant air shower properties Q . 
The analytical approach, although naturally simpli- 
fied, helped understand the important coherence ef- 
fects present because the scales in the air shower 
are comparable to the observing wavelengths of the 
radio emission. A comparison with historical data 
showed that the analytical model was also able to 
qualitatively reproduce the frequency (Fig. 1) and 
lateral (Fig. 2) distributions available at that time. 

3.2. Time-domain Monte Carlo geosynchrotron 
model (REASl) 

In a next step, the geosynchrotron model was im- 
plemented in a completely different technical way, as 
a Monte Carlo code (REASl) working in the time- 
domain [sj. Based on the same analytical parame- 
terisations for the underlying air shower properties, 
the code was able to nicely reproduce the earlier 
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Fig. 1. Electric field frequency spectrum of a vertical 10^^ eV 
air shower calculated with the analytical geosynchrotron 
model in comparison with various historical data. Solid: 
centre of illuminated area, short-dashed: 100 m from centre, 
long-dashed: 250 m. 
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Fig. 2. Radial dependence of the 55 MHz electric field com- 
ponent for a vertical 10^*^ eV air shower calculated with 
the analytical geosynchrotron model Q in comparison with 
historical data. Horizontal lines denote predicted detection 
thresholds for (from top to bottom) one antenna, an array 
of 10 antennas and an array of 100 antennas. 

findings of the analytical description. Because ex- 
perimental data to compare these results to were 
sparse at that time, this cross-check was an impor- 
tant milestone in the development of the model. In 
addition, one was now able to simulate air showers 
also for geometries other than vertical incidence and 
to produce detailed information on the polarisation 
of the radio emission. 

With the REASl code, the predicted radio emis- 
sion from cosmic ray air showers could for the first 
time be analysed in detail. Its dependence on im- 
portant air shower parameters such as shower zenith 
angle (Fig. 4), energy of the primary particle (Fig. 
3), observing frequency and magnetic field configu- 
ration as well as the polarisation characteristics of 
the radio emission (Fig. 5) were analyzed and pa- 
rameterised in a next publication [9] . In parallel, the 
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prediction that inclined air showers should have a 
more favourable footprint for radio detection than 
vertical showers had been independently made by 
Gousset, Ravel and Roy with an analytical consid- 
eration of the geometric effects arising from the air 
shower geometry pjj . 
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Fig. 3. Scaling of the REASl-simulated 10 MHz electric field 
emitted by a vertical air shower as a function of primary 
particle energy Bp Q. From top to bottom: 20 m, 100 m, 
180 m, 300 m and 500 m to the north from the shower centre. 
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Fig. 4. Contour plots of the 10 MHz absolute field strength 
for emission from a lO^'^ eV vertical air shower (left) and a 
lO"*^^ eV shower with 45° zenith angle (right) as simulated by 
REASl ti]. Contour levels are 0.25 m'^ MHz"! apart. 

When the first modern measurements of radio 
emission from cosmic ray air showers became avail- 
able at that time from the LOPES [l2| and CO- 
DALEMA iJJj] experiments, interest in the physics 
of the radio emission rose strongly, and other groups 
started their own modelling efforts. 

3.3. Frequency- domain EGS-based Monte Carlo 

Engel, Kalmykov and Konstantinov tackled the 
problem with a Monte Carlo approach in the fre- 
quency domain lj| based on a special version of 
the EGS code for the simulation of electromagnetic 
cascades. While limited by computation time to 
fairly low energies, their approach allowed a first 
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Fig. 5. Raw (unfiltered) pulses in the individual linear po- 
larization components at 200 m distance to the north-west 
from the centre of a 10^^ eV vertical air shower simulated 
with REASl p. Solid: east-west component, dashed: north- 
south component, dotted: vertical component. 



investigation of the importance of the Askaryan- 
like Cherenkov contribution, based on a realistic 
description of the refractive index profile of the 
atmosphere. Their results confirmed the expecta- 
tions that the geomagnetic emission dominates over 
the Cherenkov contributions (Figs. 6 and 7). Their 
later analyses showed that it is, however, in fact not 
straight-forward to cleanly differentiate between a 
pure geomagnetic and a pure Cherenkov compo- 
nent, because under realistic conditions the two mix 
and can no longer be clearly disentangled. 
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Fig. 6. Lateral distribution of the 30 MHz radio signal from 
a 10 TeV vertical air shower simulated with the EGS-based 
Monte Carlo code of |14fl . The geomagnetic component dom- 
inates the emission. 
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Fig. 7. Frequency spectra of 1 TeV vertical air showers 
simulated with the EGS-based Monte Carlo code of [l^ . 

3.4. Time-domain Monte Carlo geosynchrotron 
model based on AIRES 



DuVernois, Cai and Kleckner took up the concept 
of time-domain Monte Carlo simulations of geosyn- 
chrotron radiation by implementing the effect in the 
AIRES air shower simulation code The radio 
emission physics is the same as that performed in the 
REASl code. However, the AIRES-based code uses 
realistic simulated particle distributions of air show- 
ers rather than parameterisations. To make these 
simulations possible without too large a hit on com- 
putation time, thinning techniques are used, the ef- 
fect of which on the radio emission predictions are 
visible in Fig. 8. The published results [15|] were qual- 
itatively similar to those of REASl, although the 
predicted electric field strengths tended to be about 
an order of magnitude higher than those calculated 
with REASl (and later REAS2). 

3.5. Geosynchrotron spectra of finite particle tracks 

In a paper by Luo [l^ , another look at the basic 
ingredients of the analytical geosynchrotron model 
[3| was taken with an analysis of the effects of finite 
trajectories on the synchrotron spectrum of parti- 
cles defiected in the Earth's magnetic field. None of 
the currently followed approaches are, however, any 
longer based on these analytical synchrotron spec- 
tra. 
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Fig. 8. Radio emission simulated with the AIRES-based 
code of DuVernois et al. [l5^] . illustrating the influence of 
particle thinning. 

3.6. Time-domain Monte Carlo geosynchrotron 
model based on CORSIKA (REAS2) 

The REAS code made the step to realistic particle 
distributions based on individual Monte Carlo air 
shower simulations with the transition from REASl 
to REAS2 performed by Huege, Ulrich and Engel 
[l3]. In REAS2, the parameterised particle distribu- 
tions were replaced with realistic distributions simu- 
lated with CORSIKA 0. Contrary to the AIRES- 
based code by DuVernois et al., the radio emission 
routines are, however, not incorporated into COR- 
SIKA directly. Instead, the particle distributions rel- 
evant for the calculation of the radio signal are sam- 
pled with several multi-dimensional histograms over 
the course of the air shower simulation with COR- 
SIKA and saved to disk for later usage in the REAS2 
code. This approach has the disadvantage that some 
information is lost during the histogramming pro- 
cess (in particular azimuthal asymmetries in the 
air shower structure and local over- and underden- 
sities). On the other hand, the approach has sev- 
eral advantages: it allowed the authors to make a 
very gradual transition from parameterised to his- 
togrammed particle distributions, with a detailed 
analysis of how this transition changes the radio 
emission ; powerful optimisation techniques can 
be applied to cut down the needed computation time 
by more than a factor of ten; and the radio emis- 
sion simulations can be studied in detail without the 
added difficulty of shower-to-shower fiuctuations. 

A comparison of the REAS2 results with those of 
REASl showed that the overall results of REAS2 are 
similar in the frequency range of interest for radio 
detection of EAS. The most important differences 
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to REASl are a much more pronounced asymmetry 
of the air shower footprint (Fig. 9), shorter pulses 
at close distances to the air shower core, and a now 
general prediction of unipolar pulses. 

The REAS2 code can also be used to study the 
effects of air shower properties on the radio emission 
in great detail. One example is shown here in Fig. 10, 
in which the contributions of different stages of the 
air shower evolution to the radio pulses at medium 
lateral distances is illustrated. This demonstrates 
that information about the air shower evolution is 
indeed encoded in the time structure of the radio 
pulses, from which it can in principle be extracted. 
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Fig. 10. REAS2-simulated contributions of different sliower 
evolution stages to the radio pulse from a vertical lO^'^ eV 
air shower seen at 525 m north from the shower centre |l7|. 

3.7. Boosted Coulomb fields and Cherenkov 
radiation 

In a recent paper Meyer- Vernet, Lecacheux 
and Ardouin discuss radio emission from cosmic ray 
air showers described in the formalism of boosted 
Coulomb and Cherenkov fields. The authors adopt a 
simplified air shower model, and the refractive index 
is approximated as constant over the atmosphere 
with its value corresponding to that at low to mid 
altitudes. They conclude that the boosted Coulomb 
fields predict electric fields of a similar magnitude as 
those predicted by the geosynchrotron model imple- 
mented in RE AS. Other properties of the boosted 
Coulomb fields such as the exponential field strength 
decrease with lateral distance, the exponential de- 
cay of the frequency spectra and the levelling off of 
the frequency spectra for very low frequencies are in 
fact also very similar to those of the geosynchrotron 
model. Another similarity is seen in the polarisation 



characteristics illustrated with contour plots in Fig. 
11. The Cherenkov fields are expected by the au- 
thors to contribute on a smaller, but not insignifi- 
cant scale. 
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Fig. 11. Contour plots of the east-west (left) and north- 
south (right) 20 MHz electric field component of the boosted 
Coulomb field radio emission from 10* electrons and 10* 
positrons as predicted in |l9|l . 

3.8. Time-domain macroscopic geomagnetic 
emission model 

Another approach to describe radio emission from 
cosmic ray air showers has been taken by Scholten, 
Werner and Rusydi , 2Q] . They chose a more macro- 
scopic approach than the geosynchrotron particle- 
by-particle description by reviving the Kahn & 
Lerche j4| picture of transverse currents arising 
from the bulk transverse motion of electrons and 
positrons in the air shower. While they use a fairly 
simplified air shower model, they demonstrate that 
in this picture, the time structure of the radio pulses 
can be related directly to the longitudinal evolution 
of the air shower, again confirming that this infor- 
mation should be directly extractable from radio 
measurements (Fig. 12). 

One striking feature of the radio signals predicted 
by their model is that the pulses are generally bipo- 
lar. This is a consequence of the fact that their pulses 
are proportional to the time-derivative of the air 
shower profile, which is always bipolar. The bipolar- 
ity also leads to frequency spectra which decrease to- 
wards low frequencies, with the zero-frequency con- 
tribution always being exactly zero. The relation 
of these findings to the contrasting results of the 
geosynchrotron model are currently under investiga- 
tion. A reconciliation of these two very different ap- 
proaches would mark a major breakthrough in the 
understanding of radio emission from cosmic ray air 
showers. 

In a second step, Werner & Scholten have used 
their approach in combination with a more realis- 
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Fig . 9. Contour plots of the REASl-simulated (upper) and REAS2-simulated (lower) air shower emission at f = 10 MHz 
llTII . The columns (from left to right) show the total field strength, the north-south polarisation component and the east-west 
polarisation component. The vertical polarisation component (not shown here) does not contain any significant flux. Contour 
levels are 0.25 fiV m~^ MHz~^ apart, outermost contour corresponds to 0.25 /iV m~^ MHz~^. 



tic air shower evolution model, zenith angles other 
than vertical incidence and a realistic refractive in- 
dex profile to study the emission physics in more 
detail. They conclude again that the transverse cur- 
rent dominates the emission, but that there are cor- 
rections which are not negligible [2l| (Fig. f 3). They 
state that the effects of the refractive index on the 
radio emission can be strong, but so far their anal- 
ysis only applies to point sources without lateral or 
longitudinal extension. 

4. Extraction of air shower parameters from 
radio signals 

With the availability of models incorporating the 
full complexity of air shower physics, including in 
particular shower-to-shower fluctuations, it seems 
worthwhile to investigate the possibilities to extract 
air shower parameters of interest directly from radio 
signals. Aside from geometrical parameters such as 
shower core position and arrival direction, the en- 
ergy and mass of the primary cosmic ray particle 
are of particular interest. In a recent paper, Huege, 
Ulrich & Engel have used the geosynchrotron model 
as implemented in the REAS2 Monte Carlo code to 
investigate this possibility . 

They find that the lateral distribution functions of 
the (bandwidth-filtered) radio signals show a char- 
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Fig. 12. Radio pulses from a vertical 10 eV vertical air 
shower calculated with the macroscopic geomagnetic emis- 
sion model of Scholten et al. [2C| . Please note that the pulses 
at different distances have been scaled for clarity. 

acteristic behaviour. While for each shower the lat- 
eral distribution can be described with an exponen- 
tial function, the slope parameter of this lateral dis- 
tribution is related to the depth of the underlying air 
shower's maximum Xmax- More importantly, there 
exists a lateral distance region, in which all lateral 
distribution functions intersect (Fig. 14), regardless 
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Fig. 13. Radio pulse polarisation components from a 5 ■ 10^ eV air shower seen at 500 m lateral distance to the south, north 
east and west from the shower core as calculated with the model described in |2l| . 



of primary particle type and shower-to-shower fluc- 
tuations. (The location of this region only depends 
on the shower zenith angle and the radio observing 
frequency window.) This is even true for air show- 
ers with primary particles of different energies, if 
the field strength profiles are normalised with the 
energy that the corresponding air shower's electro- 
magnetic cascade deposits in the atmosphere. As a 
direct consequence, a measurement at this charac- 
teristic lateral distance would allow one to directly 
deduce the energy deposited in the atmosphere by 
the electromagnetic cascade (Fig. 15), which in turn 
can be related to the energy of the primary parti- 
cle. As shower-to-shower fluctuations do not influ- 
ence the radio signal in the intersection region, the 
intrinsic RMS spread of this energy determination 
is very low (~ 3%). It should be noted that this is a 
principle limit, not taking into account any experi- 
mental uncertainties. 

If one measures the slope of the lateral distribu- 
tion function, e.g., by combining the measurement 
in the intersection region with one at a larger lateral 
distance, one can furthermore get a handle on the 
ATniax value of the individual air shower, as shown 
in Fig. 16. The intrinsic RMS spread is around 15- 
20 g cm~^, again not taking into account experimen- 
tal uncertainties. For comparison, fluorescence mea- 
surements yield experimental resolutions of ~35- 
20gcm~^. If the value of Xmax is known for a spe- 
ciflc shower, it can in turn be related to the mass of 
the primary particle. 

While the details (in particular scales) of these 



predictions do depend on the details of the under- 
lying model, the qualitative behaviour can be con- 
sidered to be more general, as it is mostly caused by 
geometrical effects that are not critically dependent 
on the details of the emission model. 
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Fig. 14. Lateral distribution of the 32-64 MHz filtered radio 
amplitude for 10^^ eV showers coming from the south and 
observers north of the shower core |22j . Estimates of peak 
radio amplitudes that would yield a signal-to-noise ratio 
(SNR) of 10 for ideal (galactic plus atmospheric), rural and 
urban noise are marked. 



5. Conclusions 

Radio detection of cosmic ray air showers has once 
again become a very active fleld of research, and with 
it, so have simulations and theory of the radio emis- 
sion physics. By now, a number of approaches for 
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article. 
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Fig. 15. A radio measurement in the intersection region di- 
rectly yields the energy deposited in the atmosphere by the 
electromagnetic cascade on a shower-to-shower basis |22| . 
The RMS spread is only ~ 3%. 
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Fig. 16. The ratio of peak field strengths in the intersection 
region and at a larger lateral distance yields direct informa- 
tion on the max value of an individual air shower |23| . 

and technical implementations of the radio emission 
models exist. In the next few years, the models will 
have to be compared in detail with the experimen- 
tal data, which are continuously improving regard- 
ing both statistics and quality. At the same time, it 
is imperative to understand in which way the dif- 
ferent models are related to each other and to rec- 
oncile contradictions existing between different ap- 
proaches. In the meantime, the modeling efforts will 
continue to make very important contributions to 
our understanding of the radio emission physics and 
its exploitation for the study of air shower physics. 
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